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1.  INTRODUCTION 


This  report  documents  and  describes  a  validation  study  of  various  versions  of  the  Navy  Surf 
Model.  The  model  has  evolved  steadily  since  its  was  introduced  to  the  Fleet  in  1988,  but  this  is 
the  first  extensive  validation  of  the  model  using  a  set  of  comprehensive  surf  zone  measurements. 

1.1  Requirements 

Today’s  Operational  Maneuver  from  the  Sea  Strategy  requires  detailed  knowledge  of  the  complex 
littoral  battle  space.  Expeditionary  forces  must  be  prepared  to  breach  varying  layers  of  coastal 
mines  and  other  obstacles,  especially  in  the  surf  zone.  The  importance  of  surf  conditions  to  naval 
and  military  operations  are  explicitly  stated  in  several  recent  publications  by  the  Joint  Chiefs  of 
Staff(1989,  1991,  1992,  1993). 

An  essential  element  in  the  success  of  amphibious  operations  is  a  series  of  accurate  surf  forecasts 
covering  the  period  of  time  from  initial  planning  until  logistics  over  the  shore  is  completed.  The 
surf  zone  can  be  a  dangerous  place;  in  view  of  the  tremendous  number  of  people  and  equipment 
entering  this  zone  during  an  operation,  it  is  critical  that  the  on-scene  commander  know  its  natural 
hazards.  If  conditions  are  forecast  to  be  lower  than  encountered,  lives  and  equipment  may  be  put 
in  jeopardy;  forecasting  conditions  too  high  may  limit  a  commander  s  options. 

For  naval  and  military  operations,  the  surf  is  described  by  five  different,  but  related,  properties 
(Joint  Surf  Manual  [Commander,  Naval  Surface  Force,  Pacific  and  Commander  Naval  Surface 
Force,  Atlantic,  1987]):  height  of  the  highest  one-third  of  breaking  waves;  the  time  period 
between  successive  breaking  waves;  the  angle  between  breaking  waves  and  the  beach;  the  type  of 
breaking  waves;  the  longshore  (littoral)  current.  A  single,  dimensionless  number  called  the 
modified  surf  index  (MSI)  provides  a  measure  of  overall  surf  conditions.  Other  characteristics 
with  varying  influences  on  operations  are  the  surf  zone  width,  wave  length  in  the  surf  zone,  the 
average  depth  of  water  at  the  point  of  breaking,  the  direction  from  which  sea  or  swell  approaches 
the  coast,  speed  of  breakers,  the  steepness  of  breakers,  the  amount  of  spray,  variability  in  the 
height  of  breaking  waves,  and  the  expected  change  in  breaker  height  (U.S.  Naval  Oceanographic 
Office,  1969).  In  this  study,  wave  height,  longshore  current,  and  surf  zone  width  are  evaluated. 

1.2  Overview  of  the  Navy  Surf  Model 

The  navy  surf  model  software  has  been  used  extensively  throughout  the  Fleet  since  it  was 
introduced.  It  is  contained  in  the  Geophysical  Fleet  Mission  Program  Library  (GFMPL),  the 
Tactical  Environmental  Support  System  (TESS),  and  the  Mobile  Oceanographic  Support  System 
(MOSS),  and  was  developed  because  previous  surf  forecasting  techniques  were  mainly  manual 
techniques  based  on  methods  dating  to  the  1950's  which  did  not  adequately  consider  local  shallow 
water  effects.  The  model  is  the  primary  software  for  operational  surf  forecasting.  The  model  is 
also  used  for  systems  development  and  surf  climate  descriptions.  For  example,  Dowling  et  al 
(1993)  and  Salsman  et  al  (1993)  used  the  model  to  make  surf  characterizations  for  the 
development  of  the  EATD  Explosive  Neutralization  System  and  Nichols  et  al  (1996)  used  the 
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model  to  establish  engineering  requirements  for  the  Advanced  Amphibious  Assault  Vehicle 
(AAAV). 

The  surf  model  software  is  comprised  of  a  surf  zone  model  program,  a  wave  refraction  program, 
and  a  number  of  utility  programs.  In  this  study,  only  the  surf  zone  model  program  is  validated, 
and  the  term  the  surf  model  or  the  model  will  be  used  to  describe  only  the  surf  zone  model 
program.  References  to  different  versions  of  the  model  are  made  by  name,  as  described  below.  A 
description  of  the  model  and  recent  improvements  is  given  in  Section  2.0.  It  is  important  to  stress 
the  difference  between  the  surf  model  and  the  refraction  program.  The  refraction  program 
calculates  refraction  and  shoaling  of  deep  water  wave  energy  shoreward  to  a  point  outside  the 
surf  zone.  The  existing  refraction  model,  which  is  based  on  the  RCPWAVE  model  (Ebersole  et  al, 
1986),  is  now  being  replaced  by  REFDIFl  (Kirby  and  Dalrymple,  1994). 

Surf  model  inputs  are  the  following:  (1)  wind  speed  and  direction;  (2)  wave  conditions,  either  in 
the  form  of  the  wave  parameters  of  height,  period  and  direction,  or  a  complete  directional  wave 
spectrum;  (3)  a  nearshore  depth  profile  that  extends  to  a  depth  of  about  15  to  25  feet.  The  input 
wave  conditions  and  the  nearshore  depth  profile  are  the  more  important  of  the  three  inputs.  The 
depth  profile  used  in  the  model  is  computed  from  the  user-provided  depth  profile  and  tide  level. 
Model  outputs  are  the  following:  (1)  a  forecast  in  the  SURFCST  format  from  the  Joint  Surf 
Manual,  (2)  MSI;  and  (3)  an  optional  listing  of  distances  offshore,  depths,  wave  heights, 
wavelengths,  percent  breaking  waves,  and  longshore  currents  at  regularly  spaced  intervals 
throughout  the  surf  zone. 

The  following  are  versions  of  the  model  that  have  been  developed: 

1.2.1.  NAVSURF 

The  original  version  of  the  model,  here  called  NAVSURF,  is  described  extensively  in  Earle  (1988, 
1989).  The  model  has  evolved  steadily,  but  the  basic  program  structure  has  remained  unchanged. 

1.2.2.  NCSCSURF 

The  first  improvements  to  the  model  were  delivered  to  the  Naval  Oceanographic  and  Atmospheric 
Research  Laboratory  in  July  1991.  Four  improvements  were  made:  (1)  horizontal  mixing  was 
added  to  provide  better  longshore  current  forecasts  for  beaches  with  offshore  bars;  (2)  wind 
generation  of  longshore  currents  was  added;  (3)  a  better  specification  of  swell  direction  inputs 
was  added;  (4)  numerical  calculations  of  the  percentages  of  different  breaker  types  was  added. 
These  improvements  are  described  in  Earle  (1991). 

The  following  year,  the  newly  improved  model  was  further  modified  to  operate  on  a  Navy  desk 
top  computer  (DTC-2)  for  the  Naval  Coastal  Systems  Command  (NCSC),  which  was  a  Sun 
workstation  with  the  Sun  UNIX  operating  system.  NAVSURF  operates  on  a  Hewlett-Packard 
9020A  with  the  HP-UX  operating  system.  The  NCSC  version,  designated  NCSCSURF, 
incorporated  the  1991  changes  and  was  also  made  more  portable  than  NAVSURF  because  HP- 


2 


specific  features,  especially  graphics,  were  abandoned  for  simple  ASCII  text  input  and  output 
files.  The  NCSC  modifications  are  described  in  Earle  (1992). 

1.2.3.  SPE_SURF 

When  NAVSURF  was  developed,  there  was  no  capability  to  input  directional  wave  spectra  for 
automated  model  input.  Thus,  the  model  was  designed  to  internally  generate  a  directional  wave 
spectrum  by  superimposing  a  modified  Pierson-Moskowitz  wind  wave  directional  spectrum 
derived  from  input  parameters  of  wind  wave  height,  period  and  direction  and  a  narrow  banded 
swell  wave  directional  spectrum  derived  fi-om  input  parameters  of  swell  wave  height,  period,  and 
direction.  In  anticipation  of  automated  spectrum  input,  NAVSURF  internally  utilized  wave 
frequency  and  directional  bands  that  were  identical  to  those  in  the  Global  Spectral  Ocean  Wave 
Model  (GSOWM),  the  Navy’s  operational  global  wave  forecasting  model  at  the  time.  Since 
NCSCSURF  contained  this  original  limitation,  it  was  modified  by  Mettlach  et  al  (1996)  to  allow 
direct  input  of  directional  wave  spectra.  The  modified  version,  designated  SPE_SURF,  allows 
input  of  directional  wave  spectrum  files  from  the  global  and  regional  Navy  wind  driven  WAve 
Model  (WAM)  and  the  Endeco/YSI  Type  1156  Wave  Track  Buoy. 

Later  versions  of  SPE_SURF  include  options  for  the  input  of  directional  wave  spectra  collected 
fi-om  the  directional  wave  array  at  the  Coastal  Engineering  Research  Center’s  (CERC)  Waterways 
Experiment  Station  (WES)  Field  Research  Facility  (FRC)  at  Duck,  North  Carolina  and  directional 
wave  spectrum  output  files  from  the  wave  model  STWAVE  (Resio  and  Perrie,  1994). 

1.2.4  EXP  SURF 

The  first  major  overhaul  of  the  model  was  done  by  Marshall  Earle  and  Eileen  Kennelly  of  Neptune 
Sciences  Inc.  (NSI)  in  the  summer  of  1996  under  the  project  Surf  Model  Improvements  funded  by 
ONR.  Many  model  improvements  including  some  that  had  been  considered  by  NSI  in  various 
applications  of  the  model  since  1988  were  incorporated  into  a  comprehensive  revision  of  the 
software,  which  was  streamlined,  modularized,  and  made  portable  for  a  range  of  operating 
systems  and  machines.  Additionally,  the  roller  theory  of  Lippmann  et  al  (1996)  was  incorporated 
as  a  model  option.  Both  ONR  and  the  SPAWAR-funded  project  Surf  Model  Upgrade  jointly 
contributed  to  the  roller  implementation  effort. 

EXPSURF  was  an  interim  step  in  the  ongoing  development  of  SURF96.  In  this  validation  study, 
EXPSURF  was  used,  under  certain  options,  to  simulate  NCSCSURF  and  SPE_SURF;  it  was 
tested  extensively  against  standard  input  and  output  cases  from  NCSCSURF  and  SPE_SURF  and 
exact  agreement  was  found. 

1.2.5.  SURF96 

The  modularity  of  EXPSURF  provided  a  framework  for  greatly  expanding  the  range  of  options 
available  to  the  model  user  while  keeping  the  original  structure  of  the  model.  The  most  recent 
version  of  the  model,  designated  SURF96,  grew  out  of  EXPSURF  and  incorporates  several 
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theoretical  and  numerical  improvements.  A  detailed  description  of  the  model  and  its  recent 
improvements  is  given  in  Section  2.0.  There  are  situations,  including  several  that  are  not 
demonstrated  by  these  validations,  where  SURF96  provides  significantly  better  results  than 
NAVSURF.  Plans  have  been  made  to  document  STJRF96  as  an  upgrade  to  NAVSURF  while 
further  improvements  are  being  made.  Specifically,  improved  longshore  current  calculations  are 
being  examined  in  ongoing  work.  In  addition,  options  incorporated  into  SURF96  for  testing  and 
validation  are  being  evaluated  more  fully  to  select  a  standard  option  suite  for  Navy  operational 
use. 

1.3  Validation  Method 

This  study  uses  two  versions  of  the  model-EXPSURF  and  SURF96,  two  sets  of  options  for  each 
version,  and  field  measurements  from  two  surf  zone  experiments.  The  first  experiment  is  the 
National  Sediment  Transport  Study  (NSTS)  experiment  on  Leadbetter  Beach,  Santa  Barbara, 
California,  conducted  during  January  and  February  1980.  Gable  (1981),  Wu  etal  (1985),  and 
Thornton  and  Guza  (1986,  1989)  have  described  this  experiment.  The  second  experiment  is  the 
Duck  Experiment  on  Low-ffequency  and  Incident-band  Longshore  and  Across-shore 
Hydrodynamics  (DELILAH)  conducted  at  FRF  in  October  1990.  Birkemeier  et  al  (1991),  Church 
and  Thornton  (1993),  Thornton  and  Kim  (1993),  and  Smith  et  al  (1993)  have  provided 
information  about  this  experiment.  Most  of  the  directional  wave  and  tide  data  were  downloaded 
from  the  FRF  through  http://frf.wes.army.mil.  Professor  Edward  Thornton  provided  depth 
survey  data  and  processed  nearshore  wave  height  and  current  data.  The  validation  data  reduction 
methods  used  to  prepare  model  inputs  and  verification  data  are  described  in  Section  3.0. 

A  total  of  4  model  runs  were  made  using  four  NSTS  cases  and  79  DELILAH  cases.  Model 
options  and  other  details  are  explained  in  section  4.0.  Each  model  run  was  made  using 
concurrently  acquired,  available  wind  and  wave  data.  Input  depth  profiles  are  from  surveys 
conducted  on  the  same,  or  previous,  day  as  the  wind  and  wave  data.  Measurements  of  wave 
height  and  longshore  current  in  the  surf  zone  are  used  to  validate  model  output.  The  model 
outputs  were  compared  to  measurements  and  the  intercomparisons  are  presented  in  Section  5.0. 
The  accuracy  of  model-derived  surf  zone  width  and  modified  surf  index  for  19  DELILAH  cases 
are  also  presented  in  Section  5.0. 
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2.  SURF  MODEL  DESCRIPTION 


This  section  briefly  describes  the  surf  model,  including  SURF96.  The  first  subsection  gives  an 
overview  of  the  model  and  subsequent  subsections  describe  specific  model  representations  or 
calculations  with  an  emphasis  on  recent  changes  incorporated  in  SURF96. 


2.1  Overview 

The  model  is  parametric  and  one-dimensional  and  is  largely  based  on  Thornton  and  Guza  (1983, 
1986).  The  assumptions  include:  (1)  approximately  straight  and  parallel  bottom  contours  within 
the  surf  zone;  (2)  a  directional  wave  spectrum  that  is  narrow  banded  in  frequency  and  direction; 
(3)  a  Rayleigh  wave  height  distribution;  (4)  linear  wave  theory.  The  directional  wave  spectrum 
used  for  model  initialization  can  be  obtained  from  a  wave  model  or  measurements.  If  an  offshore 
directional  wave  spectrum  is  used,  it  must  be  properly  refracted,  shoaled  and  diffracted  to  a  user- 
selected  starting  depth  dstart,  which  should  lie  seaward  of  the  breaker  line.  At  dstart,  the 
directional  wave  energy  distribution  is  reduced  to  three  representative  physical  values:  (1)  the 
direction  of  the  vertically-averaged  wave  momentum  flux,  (2)  the  incident  root-mean-square 
wave  height  and  (3)  the  dominant  wave  frequency  The  direction  of  the  wave  momentum 

flux  is  computed  by  the  method  of  Higgin  et  al  (1981)  with  the  important  assumption  that  the 
directional  wave  spectrum  at  dstcivt  is  narrow  banded  in  both  direction  and  frequency.  From 
dstart,  the  model  incrementally  calculates  the  root-mean-square  wave  height  along  a  transect 
normal  to  the  beach  to  very  near  the  beach  at  the  still  water  level.  As  waves  move  through  the 
surf  zone,  the  average  rate  of  energy  dissipation  due  to  wave  breaking  and  frictional  dissipation 
balances  the  gradient  of  shoreward  energy  flux.  Energy  is  extracted  using  the  energy  dissipation 
of  a  propagating  bore  modeled  after  a  weighted  Rayleigh-distribution  of  wave  heights.  Longshore 
current  calculations  at  each  increment  are  based  on  longshore  current  theory  using  radiation  stress 
(Longuet-Higgins,  1970a,b). 

The  width  of  the  surf  zone  is  based  on  either  (1)  the  percentage  of  breaking  waves  or  (2)  the 
percentage  of  breaking  waves  and  the  location  of  maximum  energy  dissipation.  In  NCSCSURF, 
the  surf  zone  is  based  on  the  more  offshore  of  two  points:  either  the  most  seaward  point  at  which 
more  than  33  percent  of  the  waves  are  breaking  or  the  point  at  which  there  is  a  maximum  in 
energy  dissipation.  For  MSI  calculations,  the  significant  breaker  height  H^  =  ^2-  and  the 
longshore  current  Fare  the  highest  calculated  in  the  surf  zone.  The  percentage  of  each  breaker 
type-spilling,  plunging  or  surging-is  found  from  a  widely  accepted  paramaterization  of  wave 
period,  wave  height,  and  bottom  slope  and  considers  the  model-determined  breaking  wave 
probability  distribution.  Breaker  angle  is  found  from  the  direction  of  the  vertically  averaged  wave 
momentum  flux  at  dstart.  MSI  is  calculated  using  the  criteria  given  in  the  Joint  Surf  Manual. 

2.2  Wave  and  Roller  Energy 

One  of  the  important  changes  to  EXP  SURF  and  SURF96  is  the  incorporation  of  wave  rollers. 
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The  flows  of  spilling  breakers  can  be  separated  into  two  layers,  an  upper  layer  of  energy  causing 
turbulent  wave  breaking  which  rides  over  a  lower  layer  of  energy  corresponding  to  organized 
wave  motion.  The  turbulent  water  above  the  wave  is  termed  a  surface  roller.  The  original  idea  of 
such  a  two-layer  system  was  introduced  by  Longuet-Higgins  and  Turner  (1974).  The  refinements 
of  Svendsen  (1984a,b)  were  applied  by  Schaffer  et  al  (1993)  to  a  Bousinesq  model  for  waves 
breaking  in  shallow  water.  EXPSURF  and  SIJRF96  incorporate  the  model  of  Lippman  ei  al 
(1995,  1996)  which  produced  results  consistent  with  measurements  from  both  a  planar  and  a 
barred  beach. 

Calculations  of  wave  height  and  longshore  current  in  the  surf  zone  use  the  wave  energy  present  at 
a  point  offshore.  The  energy  per  unit  surface  area  in  a  wave  is  calculated  as 

(1) 


where  p  is  water  density  and  g  is  the  acceleration  due  to  gravity.  is  the  root-mean-square 
wave  height.  When  a  wave  breaks,  the  turbulent  water  spilling  off  the  front  of  the  wave  is  defined 
as  a  wave  roller.  The  energy  per  unit  area  associated  with  a  roller  is  given  as 
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where  c  is  the  phase  speed  of  the  wave, /is  the  zero  crossing  frequency,  is  the  height  of  the 
wave  at  breaking,  h  is  water  depth,  and  o  is  the  angle  the  roller  makes  with  the  body  of  the  wave 
A  value  o  equal  to  5  degrees  was  used  in  SIJRF96.  It  is  important  to  note  that  by  definition  a 
wave  roller  is  a  spilling  wave. 

2.3  Energy  Dissipation  in  the  Surf  Zone 

As  a  wave  propagates  into  the  surf  zone,  some  of  its  energy  is  dissipated.  If  wave  roller  energy  is 
not  considered,  as  in  NAVSURF,  the  energy  flux  balance  equation  is  written  as 

d{EccosQ) 


where  is  the  wave  energy,  is  the  wave  group  velocity,  0is  the  wave  direction  relative  to 
shore  normal  x,  positive  seaward,  and  H  is  root-mean-square  wave  height,  is  the  ensemble 
averaged  dissipation  function,  modeled  after  a  breaking  wave  distribution  function.  It  has  the 
form 
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(4) 


where  5  is  an  empirical  coefficient,  and  p^(H)  is  the  probability  distribution  for  breaking  waves 
described  as 

p,{H)=W{H)p{H)  (5) 


The  distribution  is  a  weighted  Rayleigh  distribution.  Weighting  functions  W(H)  have  been 
described  by  Whitford  (1992b)  and  Thornton  and  Guza  (1983).  See  Section  2.4  for  a  complete 
description  of  the  available  weighting  functions  in  SURF96. 

If  wave  roller  energy  is  considered,  the  energy  flux  balance  becomes 

djEj^^cosQ)  diEccosQ)__^^^ 
dx  dx  '' 


where  the  dissipation  term  now  becomes  a  function  of  the  roller  term  and  is  written  as 

<t>‘-egf-jCosojH^p,(H)dH 


(T) 


2A  Breaking  Wave  Probability  Distribution  Weighting  Functions 

The  solutions  of  equations  (4)  and  (7)  involve  integrating  a  weighted  Rayleigh  distribution, 

fp,(H)=fW(H)p(H)dH  (8) 


A  number  of  weighting  functions,  W(H),  have  been  described  in  the  literature.  W(H)  must  range 
from  zero  to  unity.  Five  weighting  functions  have  been  incorporated  in  SIJRF96: 


W,(H)= 


rms 

.  yh  . 


(9a) 
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where  y  is  an  empirical  value  of  0.42,  and  h  is  water  depth. 

The  user  can  select  which  weighting  function  to  use  throughout  the  model.  Note  that  in 
NCSCSURF  and  SPE_SIJRF,  equation  (9d)  was  used  as  the  weighting  function  for  the  energy 
dissipation  function,  <6b>. 

2.5  Solution  to  Wave  Height  in  the  Surf  Zone 

Equation  (6)  is  solved  using  a  numerical  forward  stepping  and  convergence  scheme  to  determine 
wave,  and  possibly  roller,  energy  along  with  values  at  each  point.  This  stepping  scheme  takes 
the  form 

(£^.CgCos0),.,  +(£  ccose),.,  =  -<e>Ax+(£^.c^cose),+(£  CCOS0).  (10) 


If  wave  rollers  are  considered,  <€>  =  <£>.  Otherwise,  <e>  =  <€^>  and  £,  =  0.0. 
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The  above  equation  is  solved  at  each  consecutive  point  stepping  through  the  surf  zone  toward  the 
beach  from  dstart.  The  right-hand-side  (RHS)  of  equation  (10)  at  dstart,  the  /th  point,  is 
determined  from  radiation  stress  parameters  calculated  from  the  initial  energy  spectra  (Earle, 
1989).  To  solve  for  the  wave  and  roller  energy  components,  and  at  the  /+!  point,  a 
convergence  routine  is  used  until  the  left-hand-side  (LHS)  of  equation  (10)  becomes  equal  to  the 
RHS.  To  do  so,  an  initial  shooting  value  of  at  point  /+1  is  set  equal  to  at  point  i  (given 
by  the  radiation  stress  calculations).  at  point  i  +1  is  iterated  until  the  LHS  equals  the  RHS 

within  a  given  tolerance.  Values  of  at  consecutive  points,  moving  toward  the  coast,  are 

calculated  in  this  way.  At  each  point,  the  incident  wave  angle  is  refracted  using  Snell’s  Law. 


2.6  Longshore  Current  Calculations 

There  are  two  options  in  SURF96  for  calculating  longshore  current  in  the  surf  zone.  The  first 
option  in  based  on  radiation  stress  longshore  current  theory  (Longuest-Higgins,  1970a,  1970b) 
using  the  balance  of  momentum  equation  in  the  longshore  direction  given  by 


^  ^  r\ 

<Xy>  +  Xy=  0 


(11) 


where  h  is  the  water  depth  and  v  is  the  longshore  current  in  the  cross  shore  x  direction.  The  first 
term  on  the  left  hand  side  is  the  horizontal  stess  in  the  along  shore  direction  exerted  by  waves  on 
the  water  and  is  given  by 


h  .  .  sin6 

X  =  <6.> - 

y  *  ^ 


(12) 


where  c  is  wave  phase  speed  and  ^is  the  angle  of  wave  apporach  with  respect  to  x  The  second 
term  represents  horizontal  mixing.  In  SIJRF96,  n  the  horizontal  eddy  viscosity  is  given  by 


/ 

p  =  2p/2^ 

\ 


p 


\  1/3 

/ 


(13) 


where  the  factor  2  is  an  empirical  constant.  The  third  term  is  the  mean  stress  due  to  bottom 
fidction  given  by 


9  \U\v 


(14) 


where  is  the  bottom  drag  coefficient  and  \U\  is  the  magnitude  of  the  near  bottom  horizontal 
wave  orbital  velocity.  Following  Thornton  and  Guza  (1986)  and  Thornton  and  Whitford  (1991) 
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\U\  is  given  by 


sinh(^//) 


(15) 


where /  is  wave  frequency,  H  is  root-mean-square  wave  height  and  k  is  the  radian  wave  number. 
The  fourth  term  is  the  wind  stress  in  the  longshore  direction  given  by 

=  Co  Pa  (16) 


where  cjs  the  wind  drag  coefficient,  ^Js  the  air  density,  \  W\  is  the  total  wind  speed,  and  W^.  is  the 
along  shore  component  of  W. 

Longshore  current  is  solved  using  a  finite  difference  approach  after  wave  dissipation  and  wave 
height  are  found  for  each  grid  point,  /  in  the  the  surf  zone  using  the  following 


K  = 


(17) 


where 


E  = 


(18) 


and 


f  A  F  ^ 


F  =- 


^x'‘ 


-  c., 


Ax' 


-{B  ,+  A  +  ,)+(^  ,  +  £  ,) 

V  /+1  i  /-!/  V  ;+!/ 


(19) 


The  terms  A„  B,  and  C,  are  given  as  follows: 


=  2p/7; 


1/3 


(20) 
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(21) 


Cjpf^  H, 

smh(kh) 


sin0 


-  \wf  sin0^ 


(22) 


where  6^  is  the  direction  of  the  wind  with  respect  to  shore  normal.  In  the  above  solution,  the 
wave  dissipation  function  depends  on  the  weighting  function  fV(H)  selected.  This  solution  for 
longshore  current  is  employed  in  SPE_SURF  and  NCSCSURF.  The  original  surf  model 
NAVSURF  does  not  contain  the  horizontal  mixing  term. 

The  second  method  for  calculating  longshore  current  employs  the  momentum  balance  equation 
without  horizontal  mixing.  Instantaneous  stress  exerted  by  the  waves  on  the  water  is  balanced  by 
bottom  stress,  which  is  a  nonlinear  quadratic  term  of  the  form 

i  =  u;u\  (23) 


where  p  is  the  density  of  water,  Cj  is  the  bottom  drag  coefficient  and  Uy  is  the  instantaneous 
longshore  current.  \U\  is  the  magnitude  of  the  instantaneous,  near  bottom,  horizontal  current 
speed,  which  is  given  by 


\IA'  =  ul  + 


[/, 


y' 


(24) 


The  components  of  U  can  be  decomposed  into  the  time  averaged  current  speed  and  the  near 
bottom  wave  orbital  speed  such  that 

U^  =  u  (25) 


and 
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In  (25)  and  (26)  the  subscript  b  denotes  near  bottom  wave  orbital  motion,  which  can  be  written 


cos  0,  =  u^.  sine 


where 


u 


H- 


+  V 


2 

b 


(27) 


(28) 


It  is  noted  that  in  this  derivation  the  mean  current  speed  in  (25)  and  (26)  is  assumed  to  be  the 
vertically  averaged  current  speed. 

Expanding  (24)  using  (27)  yields 


\U\  =  +  lU^u^cosd  +  +  2U. 


o  t  ft 


2x1/5 


(29) 


The  time  average  of  (23)  over  several  wave  periods  assuming  stationary  wave  and  current 
conditions  is  written 


<T,‘>  =  D,<  > 


<e^  sin0 
c 


(30) 


where  D2  =  /?  Cj  =  constant.  After  rearranging  (30)  to  make  the  variable  0 

Q  -  <(v  +  sin0) 


(31) 


the  equation  for  the  time-  and  vertically  averaged  longshore  current  <v>  can  be  written, 


O  -  <  D.^'^  u  sin0  > 

<v>  =  21 _ ’  _ 

<  Z),''2  > 


(32) 


Linear  wave  theory  provides  relationships  for  w„.. 


gHT 


(33) 


u 


w 


2L  cosh 


iTZh^ 

~Tj 


where 


T  = 


^  gk  tanh(^/t)^ 


An 


and  L  =  lulk  and  k  is  the  radian  wave  number. 
The  undertow  term  of  Dally  and  Dean  (1986) 


(34) 


(35) 


is  included  in  SURF96  and  has  been  tested;  however,  for  this  study,  was  set  to  zero.  Equation 
(32)  is  solved  iteratively  using  a  shooting  method  such  that  at  each  grid  point  where  H,  h  and 
are  known,  <v>  is  initially  set  to  zero,  then  the  difference  between  the  left  and  right  hand 
side  of  the  equation  is  compared  to  a  small  tolerance.  The  iterative  value  of  <v>  is  repeatedly  put 
back  in  (32)  as  Uy  until  the  tolerance  is  met.  Note  that  <U>  in  <D,  >  is  <v>. 

2.7  Determination  of  Breaker  Types:  Spilling,  Plunging,  Surging 

In  the  old  version,  the  calculation  of  the  percentage  of  spilling,  plunging  and  surging  waves  is 
based  on  a  single  point  at  the  maximum  breaker  height  location.  In  SURF96,  averaged 
percentage  values  are  used.  Inside  the  surf  zone,  where  the  bottom  slope  is  positive,  i.e.,  where 
depths  increase  with  increasing  distance  from  shore,  and  the  percent  breaking  types 
determined  fromiif,.,„^  are  ranked  from  highest  to  lowest.  For  the  highest  one-tenth  ofH^,  values, 
the  average  of  each  associated  breaker  percentages  is  determined  and  the  corresponding  values  of 
percent  of  spilling,  plunging  and  surging  waves  are  averaged.  The  averages  make  up  the  final 
output  of  the  percentage  of  breaker  types  found  in  the  surf  model  output  file. 

Earle  (1989,  1991)  describe  how  probability  distributions  of  breaker  type  are  calculated  at 
locations  in  the  surf  zone  as  a  function  of  the  probability  distribution  for  breaking  waves,  wave 
period,  and  local  bottom  slope.  Comparisons  of  analytic  and  numerical  solutions  to  equations  (4), 
(7),  and  (8)  revealed  close  agreement  at  depths  greater  than  a  few  feet;  but  in  shallower  water,  the 
two  solution  types  differed;  the  analytic  solution  sometimes  yielded  values  of  W(H)  greater  than 
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unity.  Thus,  in  SURF96,  numerical  integration  is  used. 


2.8  Self  Starting  Routine 

Since  wave  parameters  are  determined  from  the  wave  energy  at  all  points  in  the  surf  zone,  it  is 
imperative  that  energy  transformation  be  accurately  modeled,  especially  in  the  outer  edge  of  the 
surf  zone  where  few  waves  are  breaking.  When  no  waves  are  breaking,  the  roller  terms  and  the 
dissipation  in  equation  (6)  should  be  zero.  It  has  been  found  that  where  the  bottom  slope  is  small 
and  dslart  is  a  depth  where  the  percentage  of  breaking  waves  is  less  than  five  percent,  the  model 
gradually  dissipates  energy  seaward  of  the  actual  breaker  line  causing  the  waves  that  enter  the  surf 
zone  to  have  too  little  energy.  In  other  words,  starting  the  model  too  far  offshore  can  cause 
excess  energy  loss. 


A  routine  was  introduced  into  SURF96  which  prevents  excessive  energy  dissipation  outside  the 
surf  zone.  If  selected,  this  Self  Start  routine  calculates  wave  energy  and  the  percent  of  breaking 
waves  at  each  point  offshore  without  dissipating  energy  until  the  percent  of  breaking  waves 
exceeds  five  percent.  The  user  may  select  this  routine  as  an  option.  The  algorithm  for  SelfStart  is 
the  following: 


1. 

2. 

3. 


Initial  wave  parameters  at  dstart  are  determined. 

Wave  velocity  and  angle  are  calculated  at  point  /  =  2.  Wave  angle  includes 
refraction. 


Wave  energy  and  subsequently  are  calculated  at  point  /  =  2  using 

{E  c  cos0) 

{E  )  - 11 

(VOS0V, 


(36) 


4.  Given  at  point  i  =  2,  the  percent  of  breaking  waves  is  calculated.  If  it  does  not 
exceed  5  percent,  repeat  steps  2,  3,  and  4  until  the  percent  of  breaking  waves  does 
exceed  5  percent. 

5.  Once  the  percent  of  breaking  waves  to  exceed  5  percent,  continue  calculating 
wave  energies  (along  with  roller  energies  if  so  selected)  using  equation  (10). 


2.9  Directional  Energy  Spectra 

SIJRF96  allows  users  to  generate  surf  forecasts  with  seven  different  energy  spectra  types.  The 
user  can  choose  from  the  following:  (1)  internally  generated  spectra;  (2)  DELILAH  spectra;  (3) 
regional  WAM  spectra;  (4)  global  WAM  spectra;  (5)  Endeco  buoy  spectra;  (6)  REFDIFl  spectra; 


(7)  STWAVE  spectra. 

If  option  1  (internally  generated  spectra)  is  selected,  a  modified  Pierson-Moskowitz  spectrum  is 
calculated  based  on  sea  and  swell  conditions  defined  in  the  surf  model  input  file.  This  is  exactly 
analogous  to  NCSCSURF. 


2.10  Surf  Zone  Width 

Surf  zone  width  in  SURF96  is  defined  as  the  distance  offshore  to  where  10  percent  of  the  waves 
have  broken  whereas  previous  versions  of  the  model  use  the  more  offshore  of  two  points;  the 
most  seaward  point  at  which  more  than  33  percent  of  the  waves  are  breaking  or  the  point  at 
which  there  is  a  maximum  in  energy  dissipation.  In  SURF96,  if  the  percentage  of  breaking  waves 
never  exceeds  10  percent,  a  warning  message  is  displayed  to  alert  the  user.  The  change  to  a  10 
percent  threshold  is  based  on  an  investigation  jointly  funded  by  ONR  and  SPAW AR  described  in 
Section  3.2.7. 


2. 1 1  Internal  Gridding 

Prior  versions  of  the  model  calculate  depths,  distances  offshore,  and  surf  zone  parameters  at  a 
user-specified  gridding  interval.  It  has  been  found  that  wide  intervals  may  not  always  resolve 
important  features  in  the  depth  profile.  Such  aliasing  may  have  a  deleterious  effect  on  subsequent 
calculations,  especially  percent  breaker  calculations,  which  are  dependent  on  local  bottom  slope. 
To  preclude  this  potential  problem,  SURF96  uses  an  internal  grid  with  an  initial  2  feet  horizontal 
spacing.  If  the  horizontal  extent  of  the  profile  is  so  great  that  array  sizes  cannot  accommodate  a  2 
feet  spacing,  the  model  automatically  regrids  to  a  greater  interval  (3  or  more  feet).  All  calculations 
of  surf  zone  parameters  are  done  at  this  internal  interval,  but  output  is  displayed  according  to  the 
user-  specified  output  interval.  A  linear  interpolation  scheme  is  used  to  interpolate  fi-om  the 
internally  gridded  values  to  the  user  defined  output  interval. 
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3.  VALroATION  DATA  REDUCTION 


The  NSTS  and  DELILAH  data  sets  provide  input  to  the  model  as  well  as  measurements  for 
verifying  model  output.  This  section  describes  the  data  sets  used.  The  NSTS  data  set,  constituting 
four  cases  without  wind  or  directional  wave  spectra,  is  very  small  compared  to  the  extensive 
DELILAH  data  set 


3.1  NSTS  Data  Description 

The  NSTS  depth  profiles  used  were  generously  provided  by  Professor  E.B.  Thornton.  A 
comprehensive  listing  of  the  NSTS  wave  and  current  data  is  given  in  Kraus  and  Larson  (1991, 
pp.  107,  B5-B9)  and  Thornton  and  Guza  (1989,  p,  368).  All  measurements  were  68  minute 
averages. 


3.1.1.  Offshore  Waves. 

Directional  wave  spectra  were  unavailable;  therefore,  swell  wave  parameters  from  Kraus  and 
Larson  (1991)  were  used  for  wave  input.  The  depths  at  which  wave  parameters  were  collected 
are  given  in  Kraus  and  Larson  (1991). 


3.1.2.  Wind 

No  wnnd  data  were  obtained;  wind  speed  was  set  to  zero  for  all  NSTS  runs. 


3. 1. 3.  Longshore  Current  in  the  Surf  Zone. 

These  data  were  taken  directly  from  Kraus  and  Larson  (1991,  pp.  B6-B9). 

3. 1.4.  Wave  Height  in  the  Surf  Zone. 

These  data  were  taken  directly  from  Kraus  and  Larson  (1991,  pp.  B6-B9). 


3. 1.5.  Depth  Profiles 

Nearshore  depth  profiles  from  Professor  Thornton  were  used  without  modification.  Personal 
communication  with  Mary  Bristow  at  the  Naval  Postgraduate  School  confirmed  that  the  depth 
profiles  contained  tide  level 
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3.2.  DELILAH  Data  Description 


The  DELILAH  experiment  deployed  87  instruments  from  the  shoreline  out  to  the  3  meter  (10 
feet)  depth  contour,  and  collected  a  vast  amount  of  data.  The  data  collected  make  it  possible  to 
rigorously  evaluate  the  performance  of  the  surf  model.  The  following  subsections  give  details  of 
the  DELILAH  measurements  used  in  this  study. 


3.2.1.  Offshore  Directional  Wave  Spectra 


Directional  wave  spectra  used  as  input  to  the  surf  model  were  acquired  at  a  nine-element  linear 
array  of  bottom-mounted  pressure  gauges  located  on  the  8-meter  (26  feet)  contour  about  900 
meters  (0.5  nautical  miles)  offshore.  Table  1  gives  information  on  the  data  acquisition  and 
processing  of  the  directional  wave  spectra;  the  data  were  processed  by  FRF.  The  basic  analysis 
algorithm  is  the  iterative  maximum  likelihood  estimator  derived  and  described  by  Pawka  (1983). 

Table  1.  DELILAH  directional  wave  spectrum  information. 


Length  of  time  series  processed 

8192  seconds 

Data  sampling  frequency  of  time  series 

4  hertz 

Number  of  data  points  in  a  data  segment 

4096  points 

Number  of  frequency  bands  averaged 

15  bands 

Number  of  half-lapped  segments  analyzed 

10  segments 

Degrees  of  freedom  of  spectral  estimates 

160  dof 

Number  of  output  frequency  bands 

29  bands 

Width  of  an  output  frequency  band 

0.00977  hertz 

Number  of  output  direction  bins  (arcs) 

91  directional  arcs 

Width  of  an  output  direction  bin 

2  degrees 

Depth  of  measuring  array 

25.8  feet  (7.86  meters) 

3.2.2.  Wind  Measurements 

Mean  vector  wind  direction  and  mean  wind  speed  were  obtained  from  time  series  of  wind 
direction  and  wind  speed  measurements  acquired  concurrently  with  the  directional  array  data  from 
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an  anemometer  19.4  meters  (63.7  feet)  above  mean  sea  level  on  the  beach.  Since  Earle  (1991) 
states  that  the  optimum  height  of  wind  measurements  used  for  model  input  should  be  10  meters 
(32.8  feet),  a  comparison  was  made  between  the  DELILAH  wind  speed  measurements  at  19.4 
meters  and  those  same  wind  speeds  reduced  to  the  10  meter  level  using  a  logarithmic  wind 
profile,  assuming  neutral  stability.  For  the  79  cases,  a  fractional  error  of  8  percent  or  less  was 
found  between  the  19.4  meter,  and  the  derived,  10  meter,  wind  speeds.  Thus,  the  uncorrected 
19.4  meter  wind  speed  was  used  for  all  runs. 


3.2.3.  Longshore  Current  in  the  Surf  Zone 

Longshore  current  measurements  were  acquired  from  an  array  of  nine  current  meters  located 
along  a  transect  normal  to  the  shore  line  and  extending  from  near  the  water  line  to  approximately 
800  feet  (245  meters)  offshore.  Five-minute  mean  measurements  covering  the  same  time  period  as 
the  directional  array  acquisition  period  were  averaged  and  used  to  validate  model  estimates. 


3. 2. 4.  Wave  Height  in  the  Surf  Zone 

Mean  values  of  were  obtained  from  pressure  sensors  located  in  the  same  horizontal  location 
as  the  current  meters.  Five-minute  mean  measurements  covering  the  same  time  period  as  the 
directional  array  acquisition  period  were  averaged  and  used  to  validate  model  estimates. 


3.2.5.  Depth  Profiles 

Two  series  of  depth  surveys  were  used  to  construct  the  nearshore  depth  profiles  used  with  the 
model.  The  first  series  of  surveys  were  the  Deep  Sled  surveys  which  encompassed  an  area  1600 
meters  (5250  feet)  along  shore  from  near  100  meters  (328  feet)  shoreward  of  the  beach  to  1800 
meters  (5900  feet)  offshore.  Contours  from  the  8-meter  directional  array  to  the  surf  zone  sensors 
are  largely  straight  and  parallel  to  the  beach;  however,  there  is  a  deep  trough  under  the  FRF  pier. 

The  second  series  of  depth  surveys,  covering  an  area  approximately  600  meters  (1967  feet)  along 
shore  and  375  meters  (1230  feet)  offshore  to  near  the  4  meter  (13  feet)  contour  are  called  the 
Minigrid  surveys.  The  minigrid  surveys  were  made  daily  during  the  course  of  DELILAH,  however 
the  survey  of  13  October  1990  is  incomplete  because  of  high  waves  on  that  day,  thus  the  12 
October  profile  was  used  for  13  October  runs.  Figure  1  is  a  contour  plot  of  one  survey. 

The  origin  of  the  FRF  coordinate  system  is  the  intersection  of  a  shore-parallel  baseline  with  the 
southern  boundary  of  FRF  property.  Positive  directions  are  toward  340  degrees  North  along 
shore  and  toward  070  degrees  North  cross  shore.  Elevation  data  are  referenced  to  the  1929 
National  Geodetic  Vertical  Datum. 


18 


The  sensors  that  acquired  wave  height  and  longshore  current  in  this  study  are  located  along  a 


Cross— shore,  m 


DEIH-AH  romi-grid  profile  lines 


Fig^ure  2 


Figure  1  DELILAH  mini-grid  profile  lines  (from  FRF,  1991). 
cross  shore  transect  at  985  meters  along  shore,  shown  as  the  square  marks  on  Figure  1. 
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Profile  Line  Numbers 


A  single  depth  profile  cross  shore  of  the  985-meter  line  was  constructed  by  taking  all  Deep  Sled 
measurements  within  25  meters  either  side  of  the  985-meter  line.  Daily  depth  profiles  were 
constructed  by  taking  measurements  10  meters  either  side  of  the  978-meter  line  of  the  minigrid 

surveys,  which  was  closest  to  the  surf 
zone  sensors.  The  single  Deep  Sled 
line  was  merged  with  each  of  the 
minigrid  lines  by  replacing  all  Deep 
Sled  measurements  over  the  range  of 
minigrid  measurements  with  minigrid 
measurements.  Figure  2,  which  gives 
the  depth  profiles  used,  is  shown  to 
illustrate  the  dynamic  evolution  of  the 
nearshore  depth  profiles  during 
DELILAH.  The  nearshore  is  terraced 
until  9  October,  after  which  a 
pronounced  bar  is  formed.  Thornton 
et  al  (1996)  have  explored  the 

"ko  35:  4C0  450  5CC  550  600  650  700  750  800 

relgti'.'O  cross  shcro  dtstanco  (f:}  generation  of  this  bar  using  an 

TT...  j..  j  energetics-based  sediment  transport 

Figure  2  DELE.  AH  nearshore  depth  profiles  for  ten  days;  6  Oct  (large 

0);  7  Oct  (small  +);  8  Oct  (solid  line);  9  Oct  (•);  10  Oct  (large  x);  1 1  Oct  *^odel. 

(dotted  line);  12  Oct  (dash  dot  line);  14  Oct  (large  +);  15  Oct  (small  x); 

16  Oct  (small  o). 

3.2.6.  Tide  Level 

Mean  tide  level  over  the  acquisition  time  was  obtained  from  a  pressure  gauge  at  the  8-meter 
array.  Tide  level  is  used  in  the  model  to  determine  the  still  water  level.  All  versions  of  the  surf 
model  permit  tide  level  variation  for  each  model  run  so  that,  for  a  given  input  depth  profile, 
different  surf  conditions  may  be  produced  for  different  tide  levels.  A  simple  example  is  that  waves 
may  break  over  an  offshore  bar  at  low  tide,  but  may  break  closer  to  the  beach  at  high  tide. 

3.2. 7  Surf  Zone  Width 

Surf  zone  width  is  an  important  parameter  for  Navy  applications  MSI  calculations  also  indirectly 
depend  on  the  surf  zone  width,  which  must  be  determined  for  MSI  validation.  The  surf  model 
provides  a  realistic  representation  that  larger  waves  break  further  offshore  than  do  smaller  waves. 
There  is  not  a  single  location  where  all  waves  break.  Thus,  the  distance  to  the  offshore  position 
where  an  appropriate  percentage  of  waves  are  breaking  is  used  to  determine  surf  zone  width.  This 
percentage  was  refined  using  available  DELILAH  video  imagery.  Lippmann  and  Holman  (1991) 
and  Holman  et  at  (1993)  have  shown  several  applications  of  video  image  processing  to  the  study 
of  nearshore  processes  during  DELILAH.  FRF  (1991,  Appendices  B  and  C)  gives  detailed 
information  on  the  location  of  cameras  and  ground  coordinate  points.  Several  10-minute  averaged 
video  image  files  in  .jpg  format  and  software  for  displaying  and  rectifying  the  images  were 
provided  by  Professor  Rob  Holman,  Oregon  State  University  through  Dr.  Todd  Holland  of  the 
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Naval  Research  Laboratory,  Stennis  Space  Center. 

Images  acquired  at,  or  suflSciently  near,  the  same  time  as  model  run  times  were  analyzed  by 
identifying  the  locations  of  maximum  contrast  between  the  offshore  area  of  little  or  no  breaking 
waves,  represented  by  a  dark  area,  and  the  oflFshoremost  area  of  breaking  waves,  represented  by 
an  area  of  relatively  high  brightness.  The  distance  between  the  still  water  point  and  the  point  of 
maximum  contrast  is  defined  to  be  the  width  of  the  surf  zone. 


Oblique  9010120900 


Rectification  9010120900 


900  950  1000  1050 


y(m) 


Figure  3  is  a  facsimile  of 
one  the  video  images 
used.  The  upper  image 
shows  the  camera  view; 
the  box  drawn  on  the 
upper  image  is  rectified 
to  a  horizontal  plane  in 
FRF  coordinates  and 
shown  in  the  lower 
image.  The  985-meter 
line,  the  still  water  point 
at  1 15  meters,  the  point 
of  maximum  brightness 
at  194  meters,  and  the 
point  of  maximum 
contrast  at  218  meters 
are  denoted.  The  high 
brightness  either  side  of 
the  still  water  point 
shows  surf  run-up  in  the 
swash  zone.  It  is  noted 
that  the  bar  feature  as 
illustrated  in  Figure  2 
started  around  200  m 
offshore  just  inside  the 
white  area  consistent 
with  the  breaking 
image. 


Figure  3  lO-minute  averaged  video  image  from  0900  EST  12  October  1990  (top)  and 
rectified  image  (bottom).  See  text  for  details. 


21 


For  those  periods  without  adequate  video  imagery,  daily,  early  morning  visual  surf  zone 
observations  reported  by  FRF  were  used  for  surf  zone  width.  The  visual  observations  are  not 
referenced  to  a  particular  transect,  but  rather  the  distance  offshore  at  which  the  first  breaking 
waves  appear  from  the  FRF  pier.  The  video  images  are  therefore  considered  to  be  more  accurate 
validation  measurements  than  the  visual  observations. 

3.2.8  General  Wave  Conditons 

Sample  plots  of  the  DELBLAH  data  are  shown  in  Figure  4  (Birkemeier  et  al,  1991).  Wave 
condition  were  low  at  the  beginning  of  the  experiment  but  increased  after  9  October  and  then 
peaked  on  13  October.  Waves  generally  aproached  from  the  southeast  Except  on  16  October, 
most  of  the  waves  can  be  described  as  having  one  dominant  peak. 
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Figure  4.  Sample  plots  for  DELILAH  Data 
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4.  MODEL  RUNS 


Eight  sets  of  model  runs,  summarized  in  Table  2,  are  made  to  build  a  basis  for  comparison 
between  the  older  versions  of  the  model  and  recent  improvements.  EKPSURF  has  the  following 
options: 


1.  Rollers  or  not 

2.  Directional  wave  spectrum  or  wave  parameters 
SURF96  has  the  following  options: 

1.  Rollers  or  not 

2.  Directional  wave  spectrum  or  wave  parameters 

3.  One  of  five  weighting  functions 

4.  Self  start  or  not 

5  Linear  or  nonlinear  bottom  friction 

Each  option  for  EXPSURF  and  SIJRF96  was  tested,  but  only  4  combinations  from  both  versions 
of  the  model  are  used  in  this  study.  The  main  objectives  are:  (1)  to  compare  roller  with  no  roller 
runs  and  (2)  to  compare  the  effects  of  linear  and  nonlinear  bottom  friction  on  longshore  current. 
The  run  options  used  are  identified  as  follows: 

1st  letter  S  =  Wave  parameters.  (NSTS  data) 

D=  Directional  wave  spectra.  (DELILAH  data) 

2nd  letter  S  =  EXPSURF  without  rollers 
E  =  EXPSURF  with  rollers 

L  =  SURF96  with  rollers,  linear  bottom  fnction,  self  start,  and  TG83 
9  =  SURF96  with  rollers,  nonlinear  bottom  fnction,  self  start,  and  TG83. 

The  SS  and  DS  runs  provide  a  baseline  for  evaluating  model  improvements.  The  SE  and  DE 
runs  were  the  same  as  the  SS  and  DS  runs,  respectively,  except  that  the  roller  term  of  Lippmann 
et  al  (1996)  was  used  in  the  latter.  The  S9  and  D9  runs  are  the  same  as  the  SL  and  DL  runs, 
respectively,  except  that  nonlinear  bottom  friction  is  used  in  the  latter.  SURF96  has  five  different 
weighting  functions,  but  all  runs  are  made  using  the  Thornton  and  Guza  (1983)  weighting 
function  given  in  equation  (9c),  which  is  the  only  one  available  in  earlier  versions  of  the  model. 
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Table  2.  Summary  of  preliminary  validation  runs.  NSTS  denotes  data  collected  at  Leadbetter 
Beach,  Santa  Barbara  1980;  FRF  denotes  DELILAH  data;  EXP  denotes  model  EXPSURF;  SURF 
denotes  model  SURF96;  N  rum  denotes  the  number  of  runs  made;  Waves  denotes  wave  input, 
where  Tp,  dp  are  swell  parameters  and  S(f,  6)  is  a.  complete  directional  wave  spectrum;  d  start 
denotes  the  depth  at  which  model  surf  zone  calculations  were  selected  to  begin;  wind  denotes 
wind  speed  and  direction  inputs  to  the  model;  Ax  in  denotes  internal  model  grid  spacing;  Ax  out 
denotes  output  interval;  roller  denotes  use  of  roller  term;  W(H)  is  the  breaking  wave  probability 
distribution  weighting  function  used,  where  TG83  is  the  W(H)  from  Thornton  and  Guza  (1983), 
equation  (9c);  selfstart  denotes  use  of  the  self  start  option;  hjrctn  indicates  linear  or  nonlinear 
bottom  friction;  and  ID  indicates  an  identifier  for  the  set  of  runs. 


Summary  of  Runs 

Data 

NSTS 

NSTS 

NSTS 

NSTS 

FRF 

FRF 

FRF 

FRF 

Model 

EXP 

EXP 

SURF 

SURF 

EXP 

EXP 

SURF 

SURF 

N  runs 

4 

4 

4 

4 

79 

79 

79 

79 

Waves 

hoTp,0p 

VTp,e, 

S(f,0) 

s(f,e) 

S(f,0) 

S(f,0) 

dstart 

ll-13ft. 

ll-13ft 

ll-13ft 

ll-13ft 

20  ft 

20  ft 

20  ft 

20  ft 

wind 

no 

no 

no 

no 

yes 

yes 

yes 

yes 

Ax  in 

4  ft 

4ft 

4ft 

4ft 

2ft 

2ft 

2ft 

2ft 

Ax  out 

4ft 

4ft 

4ft 

4ft 

4ft 

4ft 

4ft 

4ft 

roller 

no 

yes 

yes 

yes 

no 

yes 

yes 

yes 

W(H) 

TG83 

TG83 

TG83 

TG83 

TG83 

TG83 

TG83 

TG83 

selfstart 

no 

no 

yes 

yes 

no 

no 

yes 

yes 

b_frctn 

linear 

linear 

linear 

non 

linear 

linear 

linear 

non 

ID 

SS 

SE 

SL 

S9 

DS 

DE 

DL 

D9 
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5.  RESULTS  AND  DISCUSSIONS 


The  results  of  the  model  runs  of  Table  2  are  presented  in  this  section  and  in  the  plots  in 
Appendices  A  and  B.  The  subplots  in  these  figures  are  the  depth,  percent  breaking  waves  and 
the  longshore  current  V,  respectively.  Each  figure  shows  the  results  of  two  model  runs.  The 
positions  of  the  DELILAH  pressure  gauges  are  shown  "o"  and  the  positions  of  the  DELILAH 
current  meters  are  shown  as  The  horizontal  position  of  the  sensors  are  relative  to  the  still 
water  level  The  following  sections  describe  the  results  of  the  model  in  estimating  wave  height, 
longshore  current,  surf  zone  width,  and  modified  surf  index. 

Note  in  Appendices  A  and  B  that  the  wave  height  output  from  the  SL  and  S9  runs,  and  from  the 
DL  and  D9  runs  are  identical.  This  is  expected  because  the  difference  in  bottom  friction  only 
affects  the  longshore  current  computation.  For  the  SL,  S9,  DL,  and  D9  runs,  using  Self  Start,  the 
surf  zone  width  encompasses  that  length  for  which  the  model  estimates  more  than  5  percent  of  the 
waves  to  be  breaking  or  broken.  For  the  SS,  SE,  DS,  and  DE  runs,  model  output  is  given  for  the 
entire  cross  shore  length  from  dstart  to  the  edge  of  the  water.  The  summary  scatter  diagrams  of 
model  output  are  shown  in  Appendix  C. 

5.1,  NSTS  Results 


J././Wave  Height 

The  NSTS  beach  has  a  rather  flat  beach  profile  without  any  bar  features.  From  Appendix  A,  the 
model  does  a  reasonable  job  in  wave  height  prediction  showing  the  first  shoaling  effect  and  then 
the  decay  due  to  breaking.  Columns  3  through  6  of  Table  3  show  the  accuracy  of  the  various 
model  runs  in  describing  wave  height  throughout  the  surf  zone  at  Leadbetter  Beach.  Column  3  is 
the  Pearson  product-moment  coefficient  of  correlation  between  model  wave  height  and  the 
corresponding  wave  height  measurement.  The  correlations  from  among  the  four  sets  of  runs  agree 
to  Vkithin  ,089  showing  little  significant  difference  among  model  runs. 

The  fourth  column  of  Table  3  gives  the  slope,  and  the  fifth  column  4  gives  the_y-intercept,  of  the 
linear  regression  curve  through  the  A''(jc  =  measurement,  =  model)  pairs  (See  Appendix  C).  The 
SS  runs  yield  the  best  slope  value,  and  the  SL  and  S9  runs  yield  the  best  bias  value.  The  root- 
mean-square  of  the  fractional  error  of  wave  height  is  lowest  from  the  SS  runs  and  highest  from 
the  SL  and  S9  runs. 

5. 1. 2  Longshore  Current 

The  cross  shore  profiles  of  longshore  current  measurements  contain  a  peak  between  100  feet 
offshore  and  the  water’s  edge,  but  except  for  the  3  February  case,  the  model  produces  a  current 
profile  that  increases  monotonically  from  offshore  to  the  water’s  edge.  This  problem  can  be  traced 
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to  the  formulation  of  the  longshore  current  which  employs,  for  all  options,  the  wave  dissipation 
function  <e^>.  A  detailed  investigation  of  the  complete  equations  for  <£^>  has  revealed  that  for  all 
W(H),  the  ratio  of  wave  height  to  depth  is  the  controlling  factor  for  wave  dissipation  and  thus 
longshore  current.  For  all  five  W(H)  in  SURF96  and  the  single  W(H)  in  previous  versions  of  the 
model,  it  can  be  shown  that 

HP 

<e.>  <  K—  (37) 


where  AT  is  a  constant  which  includes  all  constant  factors  in  (4)  and  (9)  and  where  p>  (}■  The 
inequality  (37)  suggests  that  in  very  shallow  water  the  rate  of  wave  dissipation  may  be  too  large, 
thus  causing  too  large  longshore  currents.  The  problem  is  unimportant  to  naval  operations  in  that 
it  only  occurs  in  depths  of  less  than  0.5  meters,  except  when  spurious  longshore  currents  at  these 
depths  are  so  large  that  that  they  mask  the  more  important,  more  accurately  modelled  currents  at 
deeper  depths.  The  problem  may  overcome  by  simply  ending  surf  zone  calculations  at  shallow 
depths.  This  deficiency  will  be  evaluated  and  eliminated  in  the  near  fixture. 

The  overall  accuracy  of  the  model  at  NSTS  was  highest  from  the  SS  (non-roller)  runs  as  shown  in 
Table  3  by  having  the  highest  correlation  between  measurement  and  model  and  the  lowest  rms 
error. 
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Table  3.  NSTS  Wave  Height  and  Longshore  Current  Verification.  Model  runs  and  data  of  3-6 
February  1980.  A/ is  the  model  version  Vis  the  number  of  model  output-measurement  pairs  for 
all  days  and  all  measurements  in  the  surf  zone  for  which  there  is  a  corresponding  model  output 
value.  R  denotes  the  correlation  coefficient  between  V  pairs,  slope  is  the  slope  of  a  best  fit  linear 
regression  through  the  N  pairs  and  bias  is  the  ^-intercept  of  the  best  fit  linear  regression  curve 
through  the  N  pairs  e  is  the  root-mean-square  of  the  fractional  error  for  N  pairs.  Subscripts  H 
denote  rms  wave  height  and  subscripts  V denote  longshore  current.  Units  of  bias„  are  feet  and 
units  of  biasy  are  knots. 


M 

N 

Rh 

slope„ 

biaSff 

Ry 

slope  y 

biaSy 

SS 

45 

0.8721 

0.9083 

0.1797 

0.2827 

0.7479 

0.6257 

0.0975 

0.4719 

SE 

45 

0.8789 

0.9055 

0.1806 

0.3042 

0.6968 

0.7844 

0.1041 

0.6008 

SL 

42 

0.8700 

0.7664 

0.1755 

0.4976 

0.5841 

0.5524 

0.1095 

0.5867 

S9 

42 

0.8700 

0.7664 

0.1755 

0.4976 

0.5270 

0.4431 

0.0954 

0.6502 

5. 2  DELILAH  Results 

At  the  beginning  of  the  experiment,  the  beach  profile  at  Duck  has  a  small  bar  feature  nearshore. 
The  increasing  wave  conditions  from  9  October  built  a  large  bar  feature  further  offshore. 

5.2.1  Wave  Height 

In  general,  the  model  gives  good  wave  height  prediction  for  both  roller  and  non-roller  runs.  The 
number  of  DELILAH  runs  makes  a  more  detailed  comparison  of  the  model  output  and 
measurements  possible  than  with  the  limited  NSTS  runs.  Table  4  gives  the  correlation  coefficient, 
slope,  and  bias  of  the  model  for  N  measurement  pairs,  but  at  each  sensor  for  the  DS,  DE,  and  DL 
and  D9  runs.  The  sensors  were  in  fixed  positions  throughout  the  course  of  DELILAH  and  due  to 
the  changing  water  level  the  model  output  nearest  the  measurement  was  found  and  used  in  the 
statistics  In  Table  4,  sensor  5  =  1  is  closest  to  the  edge  of  the  water  and  .S  =  9  is  the  most 
seaw'ard  (See  Figure  1). 

The  SURF96  runs,  i.e.  the  DL  and  D9  runs  yielded  the  highest  overall  correlation  and  a  slope  near 
unity  from  sensor  S  -2  located  in  relatively  shallow  water  where  there  was  often  a  trough 
shoreward  of  the  bar.  The  DE  runs  produced  high  correlations  and  slopes  near  unity  seaward  of 
the  bar. 

Scatter  diagrams  of  measured  and  model  wave  heights  from  the  79  DELILAH  runs  for  all 
measurements  taken  together  and  all  measurements  by  sensor  are  given  in  Appendix  C. 
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5.2.2  Longshore  Current 


Compared  to  measurements,  the  model  tends  to  underpredict  the  peak  of  the  longshore  current, 
which  is  generally  shoreward  of  the  bar,  and,  as  with  the  NSTS  runs,  to  give  inaccurate, 
monotonically  increasing  currents  close  to  the  water’s  edge.  However,  the  model  did,  at  times, 
provide  remarkably  good  current  profiles,  as,  for  example,  the  1300  EST  6  October  and  1300 
EST  7  October  cases.  Close  examination  reveals  that  all  versions  of  the  model  produce  the  best 
longshore  current  profiles  when  the  tide  is  low.  The  more  common  monotomc  increase  m 
longshore  current  near  the  water’s  edge  on  each  of  the  NSTS  runs  and  on  the  DELIL^  runs  at 
higher  tide  does  not  exist  on  the  DELILAH  runs  at  low  tide.  The  physical  reason  for  better  model 
performance  at  low  tide  appears  associated  with  the  greater  loss  of  wave  energy  when  the  top  o 
the  bar  is  very  shallow  and  thus  a  greater  dissipative  influence.  The  disscussion  of  poorly 
modelled  wave  dissipation  due  to  the  ratio  of  wave  height  to  water  depth  in  Section  5.1.2  which 
describes  NSTS  longshore  current  results,  is  also  relevant  to  the  DELILAH  runs.  From  the 
general  tendency  for  the  modelled  peak  of  longshore  current  to  be  too  low  at  all,  but  the  most 
shallow,  depths,  except  when  wave  dissipation  has  previously  occurred  at  a  shoreward,  shallow 
location,  it  is  confidently  inferred  that  wave  dissipation  is  controlled  by  more  than  the  ratio  of 
wave  height  to  water  depth. 


Table  5  summarizes  the  correlation,  slope,  and  bias  of  all  DELILAH  measurement-model  output 
pairs  by  sensor  for  all  model  runs.  The  correlation  coefficient  is  consistently  highest  for  the  DE 
runs;  however,  the  lowest  rms  errors  are  from  the  DS  runs.  Scatter  diagrams  of  the  longshore 
current  measurements  and  model  output  are  given  in  Appendix  C. 
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5,2.3  Surf  Zone  Width 


It  is  not  possible  to  validate  the  surf  zone  width  from  the  limited  NSTS  data  set,  however  both 
visual  observations  and  video  imagery  from  DELILAH  have  made  it  possible  to  validate  model 
estimates  of  surf  zone  width.  Table  6  gives  the  model  surf  zone  width  for  21  model  runs.  The 
simpler  definition  of  surf  zone  width  in  SIJRF96  produced  improved  surf  zone  estimates  in  some 
cases,  especially  those  in  which  there  were  extremely  narrow  widths  from  the  DS  and  DE  runs. 
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Table  6.  Verification  of  surf  zone  width  in  feet  from  22  DELILAH  cases.  Columns  DAYmA  TIME  denote  model  run  time  in  October 
1990  (EST)  DS  denotes  surf  zone  width  from  surf  model  runs  using  EXPSURF  simulating  SPE  SURF.  DE  denotes  runs  using 
EXPSURF  with  rollers.  DL  denotes  runs  using  SURF96  with  linear  bottom  friction.  D9  denotes  runs  using  SURF96  with  nonlinear 
bottom  friction.  Column  pairs  time-pier  and  time-video  denote  FRF  surf  zone  width  observations  (in  bold)  from  pier  and  surf  zone 
width  (in  bold)  determined  from  10  minute  averaged  video  imagery. 
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5.2.4  Modified  Surf  Index 

Although  the  MSI  cannot  be  validated,  a  consistency  check  can  be  made^  The  formula  for  MSI 
given  in  the  Joint  Surf  Manual  uses  the  maximum  wave  height  in  the  surf  zone  maximum 
longshore  current  speed  in  the  surf  zone,  breaker  angle,  wave  period,  and  breaker  type.  The 
method  used  to  check  MSI  is  to  replace  the  contribution  of  wave  height  and  longshore  current 
from  the  model  with  the  respective  contribution  from  measurement,  while  keeping  the  mode 
contributions  of  wave  angle,  wave  period,  and  breaker  type.  Wave  angle  and 
actual  directional  wave  spectra  and  there  are  no  objective  measurements  of  breaker  type^In  view 
of  the  fact  that  the  four  versions  of  the  model  for  each  case  will  give  varying  estimates  of  breaker 
type,  MSI  was  checked  for  each  version  of  the  model. 

The  19  cases  used  for  MSI  comparisons  are  given  in  Table  7,  which  lists  values  of  wave  height 
and  longshore  current  for  the  measured,  and  four  model  versions.  The  measured  values  are  the 
highest  values  in  the  surf  zone  as  determined  from  the  greatest  zone  width  for  each  run  in  Tab  e 
T^ble  7  reveals  a  general  tendency  for  the  model  to  provide  slightly  higher  wave  heights  than 
observed  and  to  provide  longshore  current  estimates  that  are  lower  than  observed. 

The  results  of  Table  7  are  used  to  compute  the  MSI  values  of  Table  8,  which  shows  that 
approximately  one-third  of  the  model  estimates  of  MSI  were  higher  than  the  mode  -measurement 
calculation  of  MSI.  Some  of  the  larger  difference  can  be  contributed  to  the  larger  longshore 
current  from  the  model.  For  instance,  the  0700  EST  8  October  case  in  Table  7  shorn  a  lonshore 
current  of  2.4  knots  as  against  the  measured  value  of  0.99.  Since  the  conmbution  of  longshore 
current  to  MSI  is  3  times  the  longshore  current  in  knots,  the  model  current  contnbuted  7.2  and 
while  measurement  contributed  2.97,  which  results  in  a  difference  in  MSI  of  4.23.  The  rms  e^or 
of  wave  height,  longshore  current  speed,  and  MSI  for  the  19  cases  are  given  in  Table  9.  Further 
discussion  is  not  useful  unless  significant  improvement  of  longshore  current  prediction  can  be 

made. 


5.3  Discussions 

Several  improvements  incorporated  into  EXPSUBF  and  SURF96  are  not  evaluated  directly  by 
these  validation  data  sets.  Although  the  Self  Start  option  is  not  explicitly  evaluated,  surf  model 
runs  at  other  beaches  show  that  starting  the  model  too  far  from  the  breaker  line  causes 
underestimation  of  wave  energy.  Also,  it  has  been  found  that  using  an  excessively  large 
calculation  interval  has,  at  times,  provided  unrealistic  results.  Inappropnate  gnddmg  is  avoided  m 
SURF96  with  the  added  feature  of  automatic  internal  gridding.  Breaker  type  has  not  been 
evaluated  ovving  to  the  lack  of  precise  measurements,  but  it  is  known  from  visual  obsemtions 
that  the  model  realistically  estimates  spilling  and  plunging  breakers.  Model  estimates  of  surging 
breakers  have  not  been  evaluated  because  such  breakers  rarely  occur.  Previous  model  versions 
calculated  the  percentage  of  breaker  types  only  at  the  point  of  maximum  occurrence,  which  made 
these  versions  excessively  sensitive  to  local  bottom  slope.  In  SimF96,  the  use  of  average  breaker 
type  over  the  entire  surf  zone  yields  estimates  that  are  less  affected  by  small  changes  in  depth. 
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The  main  objectives  of  this  validation  report  have  been  to  validate  and  test  the  robustness  of 
SIIRF96  under  a  wide  range  of  wind,  tide,  and  wave  conditions  at  both  barred  and  straight 
beaches.  In  practice,  some  of  the  cases  should  not  have  been  included  because  the  conditions  at 
the  time  violate  a  basic  assumption  in  the  model — some  DELILAH  directional  wave  spectra 
contain  two  energy  peaks,  which  violates  the  assumption  that  wave  energy  be  narrow  banded  in 
frequency  and  direction  in  the  Thornton  and  Guza  formulation.  However,  they  are  included  to 
evaluate  how  badly  the  model  will  perform  under  such  a  situation.  The  present  longshore  current 
model  is  not  supposed  to  work  well  for  barred  beaches.  But  the  results  for  barred  beaches  are 
included  for  evaluation  purpose 
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Table  8.  DELILAH  Modified  Surf  Index.  Verification  of  MSI  using  observed  surf  zone  width  and  maximum  measured  wave  height  and 
maximum  absolute  value  of  measured  longshore  current.  The  respective  contributions  of  model  wave  height  and  current  to  MSI  have 
been  replaced  by  the  respective  contribution  of  measured  wave  height  and  current.  Model  contributions  to  MSI  from  other  factors  left 
unchanged.  Capital  M  denote  best  MSI  estimates  from  model  and  measurements;  small  m  denote  model  derived  MSI.  Model  MSI  that 
is  higher  than  model-measurement  MSI  is  shown  in  bold. 
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Table  9.  Root-mean-square  of  the  fractional  error  in  the  longshore  current  and  wave  height  in  the  surf  zone,  and  MSI 
DELILAH  cases  using  four  versions  of  the  surf  model. 
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6.  CONCLUSIONS 


In  general,  the  surf  model  accurately  estimates  wave  height  throughout  the  surf  zone.  Compared 
to  older  versions,  the  roller  option  provides  better  results  for  higher  wave  conditions;  however, 
criteria  for  automatically  choosing  the  roller  option  must  be  developed  for  an  operational  version 
of  the  model  The  validation  tests  have  shown  that  the  upgraded  model  is  computationally  robust 
and  provides  reasonable  surf  zone  estimates.  The  major  deficiency  of  the  model  is  in  its  longshore 
current  computation,  which  underpredicts  the  maximum  value,  and  tends  to  give  values  that  are 
too  high  in  very  shallow  water.  This  deficiency  is  closely  associated  with  the  basic  formulation  of 
the  wave  energy'  dissipation  function.  The  implementation  of  nonlinear  linear  bottom  stress 
formulation  of  longshore  current  does  not  provide  significant  improvement  over  the  older 
formulation.  The  longshore  current  model  will  be  the  focus  for  fijture  ONR  and  SPAWAR 
projects. 

The  wide  range  of  options  in  the  model  produce  largely  similar  results  for  wave  height  and 
longshore  current;  however,  the  best  set  of  options  for  given  conditions  must  be  determined.  A 
decision  tree  approach  should  be  developed  to  assist  users  in  choosing  these  options. 

The  NSTS  and  DELILAH  data  set  are  a  valuable  resource  for  further  testing  of  the  surf  model, 
and  the  DELILAH  set  has  been  used  for  wave  refraction  model  verification.  However,  a 
deficiency  in  this  study  has  been  a  lack  of  observations  of  breaker  type.  Future  validation  data  sets 
should  include  these  important  observations  owing  their  importance  in  surf  zone  descriptions. 
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Appendix  A 

Model  and  Data  Comparison 
for  NSTS 
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Model  and  Data  Comparison 
for  DELILAH 
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Scatter  Plots 
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